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I.  INTRODUCTION 


This  report  covers  progress  in  the  development  of  the  linear  predictive 
coding/high-speed  photoconductor  sampling  switch  concept  for  wide-dynamic- 
range  digitization  of  wideband  HFDF  data,  using  a  synchronously  driven 
laser/fiber-optic  system.  The  background  of  the  concept  and  previous  work 
is  outlined  in  a  report  dated  28  March  1983  by  McKnight  and  Speiser  [1],  as 
well  as  in  patent  disclosures  under  Navy  case  numbers  66,297,  and  68,315. 
Development  activity  reported  here  is  in  the  areas  of:  (1)  simulation, 
modeling,  analysis,  and  demonstration  of  noise/jammer  cancellation  in  speech 
data  by  using  linear  predictive  coding;  (2)  assessment  of  a  finite- impulse - 
response  digital  filter  approach,  using  a  tapped  delay  line  with  fixed  tap 
weights  and  an  over-sampled  signal  to  accomplish  linear  prediction  with 
"arbitrary"  accuracy;  (3)  analysis  of  the  physical  properties  of  photo¬ 
conductor  switches  and  sample/track-and-hold  circuit  parameters  as  they 
limit  the  data  word  accuracy  in  an  A/D  converter  application;  and  (4)  cur¬ 
rent  and  future  development  activities. 


II.  APPLICATION  OF  A  LINEAR  PREDICTIVE  ALGORITHM 
FOR  REMOVING  A  JAMMER(S)  FROM  SPEED  DATA 

by  D.J.  Kaplan  and  W.H.  McKnight 


A.  THEORETICAL  BACKGROUND 

The  linear  prediction  model  is  based  on  the  use  of  an  Mth-order  linear 
predictor  of  the  sampled  (and  digitized)  signal  s(n),  where  n  denotes  the 
nth  sample,  which  requires  a  linear  combination  of  the  previous  M  samples 
(see  Markel  and  Gray  [2]).  If  s(n)  denotes  the  predicted  sample  and  e(n) 
the  prediction  error, 

M  M 

e(n)  -  ^  a^s(n-i)  -  s(n)  +  a^s(n-i)  -  s(n)  -  s(n) 

i-o  i-1 

where  „ 

M 

s(n)  -  -  ^  ais(n-i) 
i-1 

and  the  coefficients  -a^,  i  =  1,  2,  ...,  M  define  the  predictor  coefficients 

to  be  found.  The  minus  sign  is  chosen  so  that  the  error  is  based  on  a 
difference  of  two  variables,  although  this  choice  is  completely  arbitrary. 
The  total  squared  error  is  given  by 


a 


“1 

y  e2(n) 


n-n 

o 
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where  n  and  n,  define  the  index  limits  over  which  error  minimization 
o  1 

occurs.  We  can  write 


*i  r  m 


1 2  n 


1  M  M 


’2  2  ais(n'i)  '2  X  S  v(n*i)s(n*j)aj* 


n-n  i-o 
o  L 


n-n  i-o  j-o 
o  J 


and  define  c^  -  ^  s(n-i)s(n-j)  so  that 


M  M 


°  -  2  I  aicuaj  ■ 

i-0  j-O 


Minimization  of  a  is  obtained  by  setting  the  partial  derivative  of  a  with 
respect  to  a^,  k-1,2,...,  M,  to  zero  and  solving  the  set  of  M  linear  simul¬ 
taneous  equations  for  the  unknown  predictor  coefficients  {a^}. 


a  a 

I  v 


2  v 


-c  ,  ,  k-1,2 . M.  Using  the  known  parameters  c.,  ,  i  - 


0,1,...,  M,  as  defined  from  the  data,  one  can  see  that  samples  s(n)  from 
nQ-M  to  n^  are  required. 

The  two  specific  cases  most  commonly  employed  in  implementing  the 
general  linear  predictor  model  are  referred  to  as  the  covariance  method  and 
the  autocorrelation  method.  Assuming  that  a  sequence  of  N  speech  samples 

(s(n)}  -  s(o),  s(l) .  s(N-l)  is  available,  the  covariance  method  is 

defined  by  setting  nQ  -  M  and  n^  -  N-l  so  that  the  error  is  minimized  only 

over  the  interval  [M,  N-l],  and  all  N  data  samples  are  used  in  calculating 
the  covariance  matrix  elements  c^ .  In  this  case,  the  coefficient  set  { j ) 

forms  a  symmetric  semidefinite  matrix  that  will  be  singular  if  the  input 
data  sequence  (s(n)}  satisfies  a  linear  homogeneous  difference  equation  of 
order  M  or  less. 


The  autocorrelation  method  is  defined  by  setting  nQ  -  -®  and  n^  -  ®  and 

defining  s(n)  -  o  for  n  <  o  and  n  >  N.  The  coefficient  set  {c^}  then  forms 

the  elements  of  a  symmetric  positive  definite  Toeplitz  matrix,  and  the  co¬ 
efficients  can  be  expressed  in  terms  of  an  autocorrelation  sequence  as 


r(|i  -  j|) 


re 


N-l-k 


r(k)  -  s(n)s(n+k)  for  k  -  0,1,...,  N. 


n-0 

y  values  of  r(k)  for  k-0,1,...,  M  are  needed  for  the  solution. 


The  two  methods  give  identical  results  when  the  data  sequence  is 
ncated  so  that  s(n)  -  o  for  n  <  M  and  for  n  >  N  -  1  -  M.  The  two  methods 
e  similar  results  when  N  »  M.  Thus  the  essential  differences  between 
se  two  methods  of  LPC  lie  in  the  manner  in  which  they  treat  the  signal 
side  the  analysis  interval  in  the  process  of  determining  predictor  co- 
icients.  The  autocorrelation  method  requires  that  the  signal  be  set  to 
o  outside  the  analysis  interval  and  a  suitable  window,  such  as  a  Hamming 
dow,  be  used  to  reduce  the  abrupt  change  in  signal  values  occurring  at 
beginning  and  at  the  end  of  the  analysis  interval.  The  covariance 
hod  avoids  truncations  of  the  signal,  but  does  require  that  an  entire 
rix  of  covariances  be  computed  from  the  data  signal.  The  covariance 
hod  was  chosen  for  this  demonstration/simulation/analysis  because  it 
oretically  gives  exact  results  when  the  input  data  consist  of  exactly  n 
plex  sinusoids.  The  algorithm  must  be  able  to  handle  a  situation  in 
ch  there  are  jammers  present,  but  no  random  signal,  so  that  the  input 
sists  of  only  pure  sinusoids. 


CODING  AND  IMPLEMENTATION 

Even  though  the  ultimate  application  of  LPC  is  for  the  HF  band  (2-32 
),  where  sample  rates  may  be  100  megasamples  per  second  or  more,  much 
wledge  and  insight  are  to  be  gained  by  applying  the  LPC  concept  at  audio 
quencies  where  it  is  relatively  easy  to  implement.  Thus  the  specific 
1  of  this  facet  of  the  general  development  is  to  develop  the  appropriate 
ftware)  techniques (s)  and  to  show  that  a  linear  predictive  algorithm  can 
used  to  remove  a  jamming  tone  (or  tones)  from  speech  data.  This  is  some- 
t  similar  to  a  typical  HFDF  situation  in  which  the  jammer  may  be  a  strong 
entional  or  unintentional  tone  (or  tones),  which  is  (are)  coherent  over 
ny"  samples,  and  in  which  the  signal(s)  of  interest  (speech  in  this  case) 
(are)  represented  as  a  low-level  signal  that  is  relatively  incoherent  or 
erent  only  over  a  very  few  samples.  Speech  is  actually  coherent  over 
ically  about  20  ms  or  200  samples  when  sampled  at  10  kHz  (the  required 
imum  sampling  rate  given  the  speech  bandwidth  of  -5  kHz).  Speech  does, 
ever,  offer  the  advantage  of  being  intelligible  so  that  one  can  readily 
e  a  subjective  judgment  as  to  the  effectiveness  of  the  jamming  tone 
cellation  via  LPC.  In  the  case  of  A/D  conversion  for  HFDF,  the  jamming 
e  (or  tones)  would  be  retained  as  a  predicted  quantity  (also  perhaps 
luding  some  of  the  desired  signal),  whereas  in  the  present  case  of  speech 
ner  cancellation  the  predicted  jamming  signal  is  subtracted  from  the 
oming  composite  signal  and  discarded  (only  the  then -un jammed  speech  is 
ained) . 
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The  finite  coherence  of  speech  signals  requires  a  prediction  for  a 
signal  (speech  plus  jammer)  sample  that  is  far  enough  into  the  future  so  as 
to  well  exceed  the  speech  coherence  length  but  not  the  jammer  coherence 
length.  Otherwise,  the  speech  signal  would  be  accurately  predicted  by  the 
LPC  algorithm  and  canceled  (or  subtracted),  along  with  the  unwanted  jammer. 
Thus  we  need  to  employ  at  least  a  200-step  predictor.  That  is,  with  each 
new  data  sample,  our  linear  predictor  would  produce  a  data  sample  prediction 
of  at  least  200  steps  (samples)  in  the  future. 

The  procedure  for  implementing  this  LPC  concept  is  to  first  define  a 
block  of  data  from  which  to  calculate  the  filter  weights  that  best  approxi¬ 
mate  the  data.  Past  data  samples  are  thus  weighted  to  predict  a  future 
sample.  The  order,  M,  of  the  filter  denotes  the  number  of  past  samples  used 
to  predict  a  future  sample,  and  the  filter  weights  are  determined  by  mini¬ 
mizing  the  error  between  actual  data  values  and  predicted  data  values  over  a 
block  of  N  data  samples.  Of  course,  M  must  be  sufficient  for  the  algorithm 
to  predict  "several"  simultaneous  jammers.  If  the  block  of  N  samples  is  not 
characteristically  representative  of  the  data,  or  if  the  signal  (including 
noise)  statistics  change,  the  filter  weights  will  not  produce  a  very 
accurate  prediction  and  the  filter  will  need  to  be  reconfigured  (a  new  set 
of  weights  calculated).  Obviously,  it  is  not  likely  to  be  necessary  to 
reconfigure  the  filter  every  clock  cycle  (with  each  new  data  point) ,  but  the 
algorithm  for  doing  so  must  run  quickly  enough  to  keep  up  with  the  data. 
One  could  arrange  for  the  filter  weights  to  be  renewed  periodically  or 
renewed  whenever  the  prediction  error  exceeded  some  predetermined  threshold 
(such  as  when  the  dynamic  range  of  an  A/D  converter  is  exceeded).  In  the 
present  case,  the  filter  was  reconfigured  every  clock  cycle  simply  because 
it  was  easy  and  convenient  to  do  so. 

The  following  parameters  summarize  the  choices  and  operation  of  this 
LPC  application. 

1.  The  filter  order  was  chosen  as  M  -  2  for  the  case  of  a  single 
jamming  tone  and  M  =  4  for  two  jammers.  This  follows  a  general  rule  that 
requires  two  filter  orders  (degrees  of  freedom)  for  every  complex  sinusoid 
to  be  canceled. 

2.  The  value  of  N  (the  number  of  past  samples  used  in  minimizing  the 
error  for  determining  filter  weights)  was  chosen  as  8  for  the  case  of  a 
single  jammer  (M  =  2)  and  16  for  the  case  of  two  jammers  (M  -  4) .  This  also 
follows  a  general  recommendation  by  the  authors  of  the  LPC  algorithm  (Markel 
and  Gray  [2])  that  the  value  of  N  be  approximately  equal  to  4M.  Good 
results  were  obtained  by  following  this  general  guideline. 

3.  Digitized  speech  data  were  obtained  as  a  file  consisting  of  2^ 
(~64k)  samples,  produced  at  a  rate  of  10  kHz  (one  sample  every  0.1  ms).  The 
speech  consisted  of  an  adult  male  voice  with  a  British  accent  saying  "Hello 
operator,  (pause),  operator,  (pause),  hello  operator,  ..."  Although  this 
is  not  a  standard  speech  segment,  it  was  felt  to  be  sufficient  for  this 
purpose.  Signal  data  were  stored  in  an  8 -bit  (fixed-point)  video  "frame 
store  memory"  format  consisting  of  256  lines  of  data,  with  each  line  con¬ 
taining  256  data  points.  Since  we  were  simply  sampling  a  waveform,  all 
sample  values  consisted  of  real  numbers  (there  being  no  phase  or  quadrature 


information  involved) .  The  original  unprocessed  speech  data  are  shown  in 
figure  1(b).  The  pixel  values  range  from  -128  to  +127.  Sample  values  for 
line  #8  are  printed  in  floating-point  format  in  figure  1(a).  Figure  1(b) 
is  an  image  of  the  entire  speech  data  file  (all  256  lines  representing  -6.4 
seconds  of  speech) ,  in  which  the  grey  level  represents  the  amplitude  of  the 
signal.  Line  #8  occurs  near  the  beginning  of  the  first  voiced  portion  of 
the  speech  data  file. 

4.  One  or  two  digitized  (sample  rate  -  10  kHz)  jamming  tones  were 
generated  in  software  and  stored  in  a  data  file.  Figure  2  depicts  one  such 
pure  tone  at  1250  Hz,  where  the  amplitude  of  this  jammer  was  set  at  256. 
Figure  2(a)  represents  the  tone  values  for  line  #8  and  2(b)  is  the  video 
display  for  this  entire  data  file.  The  tone  amplitude  was  chosen  to  give  a 
ratio  of  jammer  to  speech  amplitude  of  approximately  2:1,  a  value  that 
results  in  the  speech  being  just  barely  audible  (although  not  intelligible) 
under  the  jammer  when  the  two  are  added  together.  This  sum  file  is  shown  in 
figure  3. 

5.  Finally,  the  sum  data  file  represented  by  figure  3  was  processed 
by  an  algorithm  written  by  Markel  and  Gray  [2],  and  the  results  are  shown  in 
figures  4  and  5.  Figure  4  shows  the  predictor  error.  This  algorithm 
contains  two  subroutines,  AUTO  and  COVAR.  As  mentioned  previously,  COVAR 
(the  so-called  covariance  method)  was  chosen  for  this  demonstration/ 
analysis.  The  driver  algorithm,  UNJAM3,  and  the  subroutine,  COVAR,  are 
shown  in  figures  6  and  7,  respectively.  UNJAM3  uses  a  second-order  (M  -  2) 
filter  and  minimizes  the  error  over  eight  samples  (N  -  8) .  The  filter  out¬ 
puts  zeros  for  the  first  eight  values  of  the  data  file,  since  initially 
there  are  insufficient  data  for  an  accurate  prediction.  As  can  be  seen  from 
figure  4(b),  the  jammer  has  been  well  removed.  When  this  file  was  played 
out  through  a  D/A  converter  and  recorded  on  audio  cassette  tape,  the  re¬ 
production  sounded  very  close  to  the  original  speech,  with  very  little 
degradation  in  intelligibility. 

Comparing  figure  4(a)  with  figure  1(a),  it  can  be  seen  that  the  values 
of  the  predictor  error  are  in  the  same  range  and  have  much  the  same  pattern 
as  the  original  speech  sample  values.  Thus,  in  this  example,  the  predicted 
values  would  not  be  outside  the  range  of  the  A/D  converter  for  the  low-order 
bits,  as  depicted  in  figures  8  and  9.  This  indicates  that  the  algorithm 
would  be  successful  in  extending  the  dynamic  range  of  A/D  conversion  in  this 
case.  The  predicted  values  shown  in  figure  5,  as  expected,  agree  well  with 
the  jammer  tone  values  shown  in  figure  2(a). 
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Figure  1 .  Original  digitized  speech  data. 
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Figure  2. 1250-Hz  jamming  tone. 
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Figure  3.  Speech  data  plus  jamming  tone. 
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Figure  5.  Predicted  data  (jamming  tone). 


LINE  #8 


PREDICTED  VALUE 


v  f. 


1 v V *  y 1 


’  $*•  <" 


C  This  prograe  UNJAfls  or  reeoves  an  interfering 

C  jaeeing  tone  fro*  a  digitized  speech  file. 

PROGRAM  UNJAM3 

DIMENSION  I(8),A1 (256),A(3),6RCI2),E(256),P(256) 
OPEN ( UNI T»1 0 , ACCESS= ' DIRECT  * , RECORDS 1 2E=2S6, 

1  MAIREC*256, TYPE*’ OLD’) 

OPEN  IUNIT*1 1,  ACCESS*’  DIRECT’,  RECDRDSIZE*256, 

1  nAIREC*256,TYPE*’NEN’> 

QPEN1UNIT*12, ACCESS*’ DIRECT’, RECCRDS1ZE*256, 

1  MAIREG*256,TYP£*’N£N’ ) 

M=2 

N*8 

C  output  zeroes  in  the  first  8  positions 
DO  1  J=1 , 8 
E(J)=0.0 
P(J)=0.0 

1  CONTINUE 

C  coepile  the  first  8  saaple  points 

READdO’DAl 
DO  2  J*i,B 
I(J)=A1(J) 

2  CONTINUE 

C  do  the  first  line 

1*1 

C  update  the  saiple  vector 

DO  4  J=9,256 
DO  3  K*l,7 
I(K)=I(K+1) 

3  CONTINUE 
I(8)=A1(J) 

C  calculate  the  filter  coefficients 

CALL  COVAR (N, I , M, A, ALPHA, SRC) 

C  predict  the  next  sasple  and  calculate 

C  the  predictor  error. 

C  PRINT  10,  J 

C  PRINT  9,A(2),AI3) 

C  PRINT  9, 1 17), 1 16) 

P(J)*-A(2)m7)>A<3)«(6) 

E(J)*I(8)-P(J) 

C  PP=P(J) 

C  EE*E(J) 

C  TYPE  10, J 

CIO  FORMAT  flN,’J*’,I12) 

C  PRINT  9, PP, EE 

C9  F0RNATdI,2E16.6> 

4  CONTINUE 
NRITElll’DP 
WRITE ( 12’ I) E 
TYPE  5,1 

5  F0RBAT<’*’,U2) 

C  80  TO  6 


C  begin  tain  loop 

DO  6  1*2,256 
READdO’DAl 
DO  7  J=l,256 
DO  8  K*l,7 
I(K)*I (K+l) 

8  CONTINUE 

1(8)*A1(J) 

C  calculate  the  filter  coefficients 

CALL  COVAR (N, I, N, A, ALPHA, GRC) 

C  predict  the  next  saepie  and  calculate 

C  the  predictor  error. 

P(J)=-A(2)tI(7)-A(3)l)d6) 

E  (J)=X (8>-P 1 J) 

7  CONTINUE 

NRITEdl’DP 
NRITEd2’I)E 
TYPE  5,1 
6  CONTINUE 

END 


Figure  6.  Driver  program. 
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SUBROUTINE:  COVAR 

A  SUBROUTINE  FOR  IMPLEMENTING  THE  COVARIANCE 
METHOD  OF  LINEAR  PREDICTION  ANALYSIS 


SUBROUTINE  COVARiN,  I,  M,  A,  ALPHA,  GRC) 

INPUTS:  N  -  NQ.  OF  DATA  POINTS 

X  <N)  •  INPUT  DATA  SESUENCE 
M  -  ORDER  OF  FILTER  <N<21,  SEE  NOTED 
OUTPUTS:  A  -  FILTER  COEFFICIENTS 
ALPHA  -  RESIDUAL  •ENERGY* 

A  -  FILTER  COEFFICIENTS 

GRC  -  'GENERALIZED  REFLECTION  COEFFICIENTS', 

IFROSRAM  LIMITED  TO  M=20,  BECAUSE  OF  THE  DIMENSIONS 
B(MI(Hd)/2),  DETAIN),  AND  CC(Md) 

DIMENSION  1(1),  Ad),  GRC d ) 

DIMENSION  6(190),  BETA (20) ,  CCI21) 

MP  *  M  t  1 
MT  *  MPJM/2 
NT  =  (HPIH)/2 
DO  10  J*1,HT 
B(J)  *  0. 

10  CONTINUE 
ALPHA  =  0. 

CCd)  =  0. 

CCI2)  -  0. 

DO  20  NP=MP,N 
NP2  =  NP  -  1 

ALPHA  =  ALPHA  ♦  KNP)tX(NP) 

CCd)  =  CCd)  +  I(NP)IKNPl) 

CC 12)  =  CC(2)  *  KNPmi(NPl) 

20  CONTINUE 
Bd)  =  1. 

BETAd)  «  CC  (2) 

GRCfli  =  -CC(i;/CC(2) 

Ad)  *  1. 

A  (2)  =  GRCd) 

ALPHA  =  ALPHA  t  SRC  d)  ICC  (1) 

MF  =  N 

DO  130  MINC*2,MF 
DO  30  J=i,HINC 
JP  *  MINC  ♦  2  -  J 
N1  =  HP  ♦  2  -  JP 
N2  s  N  +  1  -  MINC 
N3  *  N  ♦  2  -  JP 
N4  =  HP  -  MINC 

CC(JP)  «  CC(JP-l)  ♦  KNDIKNl)  -  I(N2)II(N3) 

30  CONTINUE 
CCd)  *  0. 

DO  40  NP*HP,N 
Mil  =  NP  -  MINC 
CCd)  *  CCd)  ♦  KNl)IX(NP) 


40  CONTINUE 

MSUB  «  (MINCIMINC-HINCl/2 
MM1  *  MINC  -  1 
N1  -  MSUB  ♦  MINC 
BIND  -  1. 

DO  80  IP*1,HM1 
I SUB  *  (IPIIP*IP)/2 
IF  (BETA! IP) )  ISO,  ISO,  SO 
SO  GAM  -  0. 

DO  60  J*1,IP 
N1  *  I SUB  ♦  J  , 

GAM  «  SAM  ♦  CC(J+1)IB(N1) 

60  CONTINUE 

GAM  =  SAM/BETA  dP) 

DO  70  JP=1,IP 
HI  «  MSUB  ♦  JP 
N2  *  I SUB  *  JP 
B(N1)  *  BIND  -  6AM4BIN2) 

70  CONTINUE 

BO  CONTINUE 

BETA (MINC)  =  0. 

DO  90  J=1,M!NC 
N1  =  MSUB  +  J 

BETA(MINC)  «  BETA  (MINC)  ♦  CCIJd)IB(Nl) 
90  CONTINUE 

IF  (BETA (NINO)  ISO,  ISO,  100 
100  S  =  0. 

do  no  ip=i,nm 

S  *  S  ♦  CC(IP)IAdP) 

119  CONTINUE 

GRC (MINC)  =  -S/BETA (MINC) 

DO  120  IP*2,MINC 
M2  *  MSUB  ♦  IP  -  1 
AilP)  *  A(IP)  ♦  GRC (MINC) IB (M2) 

120  CONTINUE 

A(MINCd)  «  SRC  (MINC) 

S  =  GRC (MINC) IGRC (MINC) IBETA (MINC) 

ALPHA  •  ALPHA  -  S 
IF  (ALPHA)  ISO,  150,  130 
230  CONTINUE 
140  RETURN 
150  CONTINUE 


C  WARNING  -  SINGULAR  MATRII 
C 

C  PRINT  9999 
C9999  FORMAT  (34H  WARNING  - 
SO  TO  140 
END 


SINGULAR  MATRIX  -  COVAR) 


Figure  7.  Covariance  method  subroutine. 


C.  RESULTS  AMD  CONCLUSIONS 


During  this  LPC  voice  implementation  study/analysis,  the  following 
parameters  were  varied.  In  all  cases,  the  jamming  spectrum  was  constant 
during  the  entire  data  file. 

1.  M  -  the  order  of  the  filter  (the  number  of  past  data  samples  used 
in  the  prediction  calculation) . 

2.  N  -  the  number  of  data  samples  over  which  the  error  is  minimized. 

3.  The  number  of  jamming  tones-1  or  2. 

4.  The  frequency  of  the  jammer(s). 

5.  The  amplitude  of  the  jammer(s)  relative  to  that  of  the  speech. 

6.  P  -  the  number  of  steps  ahead  predicted  by  the  algorithm  (single 
step  or  multistep). 

7.  The  frequency  with  which  the  filter  coefficients  are  recalculated. 

Results  related  to  each  of  the  above  factors  were  as  follows: 

1.  As  is  known  theoretically,  M  must  be  at  least  twice  the  number  of 
complex  sinusoids  in  the  input  data.  This  was  observed  experi¬ 
mentally.  When  there  were  two  jamming  tones  (one  at  1250  Hz  and 
one  at  1406  Hz)  added  to  the  speech  data  and  M  was  set  at  2,  the 
algorithm  was  unable  to  remove  both  tones  effectively  and  a 
buzzing  tone  was  heard  throughout  the  file  corresponding  to  the 
beat  frequency  between  the  two  jammers.  With  M  -  4,  the  buzzing 
was  eliminated  but  the  intelligibility  of  the  speech  was  less  than 
for  the  M  -  2  case.  This  result  occurs  because  speech  can  be  well 
approximated  as  a  sum  of  M'  pure  sinusoids,  where  M'  is  typically 
on  the  order  of  10-20.  Thus,  as  M  is  increased,  the  LPC  algorithm 
becomes  a  better  predictor,  not  only  of  the  jammer(s)  but  also  of 
the  speech  as  well.  The  predictor  error  therefore  becomes 
smaller.  This  is  a  consequence  of  speech  not  being  well  modeled 
in  this  case  as  a  random  residual  since  it  has  a  finite  coherence 
length.  However,  for  the  HFDF  application,  the  residual  should  be 
well  modeled  as  a  more  truly  random  variable,  and,  in  that  case, 
the  order  of  the  filter  should  be  large  enough  to  handle  the 
number  of  degrees  of  freedom  that  the  jammers  represent. 

2.  As  stated  by  Markel  and  Gray  [2],  N  must  be  at  least  four  times 
greater  than  M  for  good  results .  This  was  observed  to  be  true  in 
our  results.  However,  N  must  be  kept  as  small  as  possible,  con¬ 
sistent  with  this  requirement,  to  keep  calculation  times  to  a 
reasonable  level.  It  was  found  that  it  took  8.5  minutes  to  run 
the  LPC  program  on  a  64 -k  sample  data  file  with  M  -  2  and  N  -  8 
(single  jammer)  when  the  filter  coefficients  were  recalculated 
every  clock  cycle.  If  one  does  not  reconfigure  the  filter  every 
clock  cycle  (and  there  is  no  need  to  if  the  jammer  spectrum 
remains  reasonably  constant),  one  can  increase  both  M  and  N  with- 
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out  adversely  affecting  the  processing  time.  Theoretically,  M 
should  be  made  large  enough  that  M  samples  will  represent  a  time 
span  that  is  large  compared  to  the  coherence  time  of  the  residual. 

The  effect  of  increasing  the  number  of  jamming  tones  was  to 
require  increasing  M,  as  explained  above.  If  the  jammers  are 
bands  rather  than  lines,  additional  requirements  will  be  placed  on 
the  LPC  algorithm.  This  case  has  not  yet  been  investigated. 

Changing  the  frequency  of  the  jammer  had  some  effect  on  the 
intelligibility  of  the  unjammed  speech.  When  the  jammer 
frequency,  f,  was  set  to  664  Hz,  the  intelligibility  of  the  speech 
was  degraded  because  f  was  close  to  the  fundamental  frequency  of 
the  speech  (-500  Hz).  This  occurs  mainly  because  a  significant 
part  of  the  speech  spectrum  is  falling  into  the  passband  of  the 
filter  and  hence  is  being  subtracted  along  with  the  jammer. 
Again,  changing  the  frequency  of  the  jammer  should  not  present  a 
problem  for  truly  random  or  broadband  signals. 

Varying  the  amplitude  of  the  jammer (s)  relative  to  that  of  the 
speech  did  not  seem  to  have  much  effect  on  the  intelligibility  of 
the  unjammed  speech  as  long  as  that  amplitude  was  no  lower  than 
about  2:1.  This  was  to  be  expected  because  the  LPC  algorithm 
works  on  the  power  in  the  coherent  part  of  the  signal.  Since  the 
speech  is  also  highly  coherent  over  the  filter  length,  it  is  not 
surprising  that  the  LPC  algorithm  would  break  down  (partially 
predict  the  speech  signal  and  subtract  it  along  with  the  jammer) 
for  low-amplitude  jammers.  Jammer- to -speech  amplitude  ratios  as 
large  as  256:1  were  tried  and,  with  proper  values  of  M  and  N,  it 
was  possible  to  accurately  reproduce  the  speech  to  an  intelligible 
degree . 

P  was  varied  in  an  attempt  to  overcome  the  problems  associated 
with  the  coherence  of  the  speech.  Values  used  were  P  —  1,  32,  and 
2000.  At  P  -  2000,  the  speech  at  2000  samples  ahead  should  have 
been  completely  decorrelated  from  its  present  value.  However,  not 
much  improvement  was  noticed  in  the  intelligibility  of  the  speech 
and,  in  fact,  an  additional  problem  was  observed  since  an  echo  of 
the  jammer  was  produced  at  the  end  boundaries  of  speech  segments. 
Thus  multistep  prediction,  at  least  for  large  values  of  P,  does 
not  seem  to  be  appropriate  for  this  application.  It  would  be  even 
less  appropriate  for  the  HFDF  application,  in  which  the  coherence 
time  of  expected  residuals  is  very  short  and  the  high  sample  rate 
will  not  allow  the  slow  algorithms  required  for  multistep 
prediction . 

The  times  at  which  the  filter  is  to  be  reconfigured  were  found  to 
be  very  important  in  the  problem  of  removing  a  jammer  from  speech. 
If  this  type  of  LPC  technique  were  to  be  used  for  HFDF  signals,  it 
would  also  be  important.  If  the  filter  is  reconfigured  during  a 
period  when  the  signal  statistics  are  nonstationary,  as  during  a 
boundary  between  speech  and  silent  segments,  the  filter  will  be 
reconfigured  incorrectly  and  the  prediction  error  will  exceed  the 


dynamic  range  of  Che  residual  A/D  converCer.  This  was  born  out  by 
observation  in  this  situation. 

In  conclusion,  the  following  observations  become  apparent  for  applica¬ 
tion  of  this  L PC  technique  to  HFDF  data: 

1.  The  filter  order  (M)  and  the  minimization  length  (N  samples) 

should  be  kept  as  short  as  possible  to  accomplish  the  objectives 

of  digitizing  the  total  signal  to  the  required  accuracy. 

2.  Single-step  prediction  should  be  used  (if  possible)  for  greater 
prediction  accuracy. 

3.  The  filter  weights  should  not  be  recalculated  every  clock  cycle 

but,  rather,  only  when  the  residual  exceeds  the  range  of  the  cor¬ 

responding  A/D  converter. 


III.  REAL-TIME  SIGNAL  PREDICTION  FOR  HF  A/D  CONVERSION 

by  J .M.  Speiser 


A.  LINEAR  PREDICTION  TECHNOLOGY 

We  have  previously  addressed  the  use  of  linear  prediction  to  increase 
the  dynamic  range  and  precision  of  analog- to -digital  converters  [1].  It  was 
shown  that  the  dynamic  range  or  precision  could  be  improved  (in  units  of 
power)  by  a  factor  equal  to  the  ratio  of  the  process  variance  to  its  one- 
sample  prediction  error  variance. 

Two  principal  concerns  in  applying  this  technique  are  (a)  ideal 
predictor  performance,  and  (b)  real-time  implementation  of  the  predictor 
when  the  sampling  rate  is  high. 

For  a  wide-sense  stationary  random  process  with  spectral  density  func¬ 
tion  s(f),  the  dynamic  range  improvement  or  ratio  of  input  random  process 
variance  to  prediction  error  variance  is  the  ratio  of  the  arithmetic  mean  of 
s(f)  to  its  geometric  mean  [3]. 

Initial  candidates  for  the  predictor  realization  were  as  follows. 

1.  Recursively  updating  an  estimate  of  the  sample  covariance  function 
and  solving  the  normal  equations  by  using  Durbin's  algorithm  [4].  This 
would  require  on  the  order  of  N- squared  multiply-adds  per  prediction  when  N 
past  samples  are  used— not  feasible  for  high-speed  real-time  implementation. 

2.  Nonreal-time  estimation  of  the  covariance  function  and  solution  of 
the  equation  for  prediction  coefficients.  This  requires  the  assumption  of 
stationarity  over  an  extended  period,  and  could  therefore  result  in  un¬ 
satisfactory  tracking  if  the  stationarity  period  was  less  than  the  sum  of 
the  covariance  update  time  and  the  prediction  equation  solution  time. 


3.  Use  of  an  adaptive  lattice  filter  (preferably  a  least-squares 
lattice  rather  than  a  gradient  descent  lattice  filter)  to  provide  order  of  N 
parallelism  in  solving  the  prediction  problem  [5],  However,  the  least- 
squares  lattice  is  a  relatively  complicated  structure,  and  requires  addi¬ 
tional  computations  to  update  the  reflection  coefficients. 

In  summary,  the  previously  proposed  predictors  have  two  areas  of  diffi¬ 
culty:  (a)  complexity  and  consequent  difficulty  of  providing  a  real-time 

implementation  with  a  latency  of  less  than  the  sampling  interval,  when  the 
sampling  rate  is  high;  and  (b)  difficulty  of  guaranteeing  prediction  per¬ 
formance  without  prior  knowledge  of  the  spectrum  of  the  input  process . 

It  appears  that  both  of  these  difficulties  can  be  circumvented  by  using 
some  little -known  results  in  prediction  theory  attributable  to  Brown  [6]  and 
Splettstosser  [7]. 

Brown  showed  that  any  deterministic  signal  or  wide -sense  stationary 
random  process  sampled  at  a  rate  above  the  Nyquist  rate  can  be  approximated 
arbitrarily  well  by  a  linear  combination  of  its  past  samples  (see  his  ref 
5).  In  fact,  he  provides  two  choices  of  such  weights  that  may  be  computed 
without  knowing  the  signal  spectrum.  Only  a  knowledge  of  the  oversampling 
ratio  is  required.  One  set  of  weights  is  given  explicitly  for  the  case  when 
the  sampling  is  at  a  rate  exceeding  twice  Nyquist: 

n 

Xn  '  1  *<t-kT> 

k=l 

where 

k+1 

akn  =  <'1)  (cos  *T) [n!/k! (n-k) ! ] , 

T  is  the  normalized  sampling  interval,  and  n  is  the  number  of  past  samples 
used  in  the  prediction.  For  this  method,  it  is  required  that  T  <  1/2. 

Brown  gives  a  second  method  for  choosing  the  prediction  coefficients 
for  any  sampling  rate  above  Nyquist  (i.e.,  T  <  1)  as  the  solution  of  a  set 
of  linear  equations  requiring  only  knowledge  of  the  sampling  interval  and 
number  of  past  samples  to  be  used  in  the  predictor: 

n 

2  akn  q(k-p)  -  q(p)  for  p  -  1, 2 , . . . ,n 

k-1 

where 

q(p)  -  sin(»rpT)/(?rpT)  . 

Brown  shows  that  for  either  of  the  above  sets  of  prediction  coeffi¬ 
cients,  the  transfer  function  of  the  prediction  error  filter  converges 
uniformly  (in  frequency)  to  zero  as  the  number  of  past  samples,  n,  becomes 
large.  This  guarantees  that  the  mean-squared  prediction  error  converges  to 
zero.  For  the  second  set  of  predictor  coefficients,  he  also  shows  that  the 


integral  of  the  magnitude  squared  of  the  prediction  error  filter  transfer 
function  is  minimized  for  each  choice  of  the  number  of  past  samples 
utilized,  thus  minimizing  the  Cauchy -Schwarz  bound  for  the  mean- squared 
prediction  error. 

Splettstosser  exhibits  a  closed- form  solution  for  asymptotically  good 
predictor  weights  for  any  sampling  rate  above  1.5  Nyquist  (i.e.,  T  <  2/3), 
but  does  not  show  a  desirable  extremal  property  for  his  weights  for  each 
value  of  n,  the  number  of  past  samples  used. 

The  reason  that  Brown  and  Splettstosser  are  able  to  obtain  perfect  pre¬ 
diction  via  oversampling  is  that  the  spectrum  of  an  oversampled  signal  has  a 
region  of  nonzero  extent  where  it  is  zero.  The  prediction  error  results  in 
Hannan's  book  [2]  then  tell  us  that  perfect  prediction  is  possible.  More 
specifically,  when  the  (un-normalized)  sampling  interval  T  is  less  than 
1/(2W),  the  spectrum  is  zero  between  W  and  the  folding  frequency  1/(2T). 
Since  the  signal  spectrum  is  lowpass  with  cutoff  frequency  W,  and  if  the 
prediction  error  filter  is  chosen  to  be  a  highpass  filter  with  cutoff 
frequency  between  W  and  the  folding  frequency  1/(2T),  the  output  of  the  pre¬ 
diction  filter  has  zero  error.  The  closer  1/(2T)  is  to  W,  the  sharper  the 
cutoff  characteristic  required  for  the  prediction  error  filter — and  hence 
the  larger  the  number  of  taps  that  will  be  required. 

The  above  discussion  shows  that  when  the  signal  is  oversampled,  there 
is  a  great  deal  of  flexibility  in  the  choice  of  prediction  filter  with  small 
prediction  error.  For  application  to  the  predictive  A/D  converter,  it  is 
desirable  that  the  predictor  be  robust  with  respect  to  perturbation  in  its 
input  data,  since  the  prediction  process  will  have  to  start  up  by  using  only 
coarsely  quantized  past  samples.  In  view  of  the  difficulty  of  analyzing  a 
feedback  loop  with  imbedded  nonlinearity,  it  will  be  necessary  to  simulate 
different  choices  of  prediction  filter  coefficients. 

Since  it  is  crucial  for  the  predictive  A/D  application  to  keep  the 
latency  less  than  the  sampling  interval,  a  strong  candidate  implementation 
is  a  hybrid  digital/analog  transversal  filter  employing  digital  multi¬ 
plication  but  analog  summation.  While  digital  summation  is  preferable  with 
respect  to  accuracy,  a  binary  tree  of  adders  to  sum  n  terms  has  a  latency  of 
log^n  addition  times.  Even  though  the  throughput  is  adequate,  the  latency 

is  inadequate  when  the  sampling  period  is  comparable  to  an  arithmetic  opera¬ 
tion  time. 


B .  MATHEMATICS 

1.  Hannan's  prediction  error  summary  used  the  spectrum  specified  in 
terms  of  radian  frequency.  It  is  both  simpler  and  more  convenient  to  use 
those  results  converted  to  ordinary  (normalized)  frequency. 

Let  the  correlation  function  for  a  random  sequence  be 
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-1/2 


i2irft  ...  .. 

e  s(f)  df 


where  s(f)  is  the  spectral  density  function  of  frequency  f. 


estimation  situation 

no  knowledge 


variance  of  prediction  error 
1/2 

J  s(f)  df  [arithmetic  mean  of  s(f)] 

-1/2 


knowledge  of  all  past 


exp 


’  1/2 

|  log[s(f)]df 

.-1/2 


[geometric  mean  of  s(f)] 


knowledge  of  past 
and  future 


1/2 


I  s(f)df 


-1/2 


[harmonic  mean  of  s(f)' 


2.  Oversampling  (deterministic  case): 

CO 

Let  h(t)  -  ^  S(t-mT) 
m--« 


H(f  ) 


|  h(t)  e'12*ft  dt  -  (1/T)  ^  «[f-(n/T)] 


n--® 


(See  M.J.  Lighthill,  Fourier  Analysis  and  Generalized  Functions,  Cambridge 
University  Press,  1964,  pp.  67-68) 


Let  g(t) 


-I 


G(f )  ei2*ft  df 


W 


g 


a(t)  -  J  G(f )  e12,rft  df 


-W 


f  _  i2»rft 

-  I  Gj(f)  e  df 


where 
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G^f) 


G(f),  |f|  <  W 
0  ,  |f|  >  V 


and  g^(t)  is  a  lowpass  filtered  version  of  g(t). 

Let  g?(t)  -  g-i(t)  h(t)  -  sampled  version  of  g,  (t)  with  sampling 

interval  T. 

00 

G2(f)  -  G1(f)  *  H(f )  -  (1/T)  ^  G^f-Cn/T)]. 

n--oo 

The  situation  is  illustrated  pictorially  below,  assuming  T  <  (1/2W)  so 
that  g2(t)  is  oversampled. 


G2(f)  is  periodic  with  period  1/T.  The  frequency  1/(2T)  is  called  the 
folding  frequency. 


f(l/T)G  .(f),  | f |  <  W 

,<f)  -  I 

[  o  w  <  I  f I  <  <: 


L  0  W  <  | f |  <  (1/2T)  . 

Note  the  following  spectral  density  function  for  an  oversampled,  wide- 
sense  stationary  (wss)  random  process.  Let  x(t)  be  a  wss  random  process 
with  autocorrelation  function  r(t)  and  spectral  density  function  s(f). 


r(t)  -  E[x(u)  x(u+t) ] 


-  J 


+i2wft 

s(f)  e  df 


where  E  denotes  statistical  expectation. 
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Let  x^(t)  be  a  lowpass  filtered  version  of  x(t),  bandlimited  to 

with  corresponding  autocorrelation  function  r..(t)  and  spectral  density  func¬ 
tion  s^(f ) . 


^(t)  -  E^x^u)  xx(u+t)|  -  |  s(f)  ei2,rft:  df  -  J 


i2*ft 

If)  e  df 


f  s(f),  | f |  <  W 

*,(f)  “  { 

[  0  ,  |  f  |  >  W 


Next,  consider  the  random  sequence  obtained  by  sampling  x..(t)  at  multi¬ 
ples  of  T,  where  T  <  1/(2W) . 

W 

E  |Xl(kT)  Xl[(k+n)T]j  -  r^nT)  -  J  s(f)e12*fnt  df  . 

-W 

Since  T  <  1/(2W),  it  is  also  true  that  W  <  1/(2T) . 


Define  82(f) 


f  s(f),  | f |  <  W 
1  0  ,  W  <  |  f  |  < 


1/(2T) 


1/(2T) 


,  f  /c\  i2irfnT 

:x(nt)  -  I  s2(f)  e  df  . 


-1/(2T) 

In  other  words,  the  spectral  density  function  for  the  oversampled 
random  process  is  zero  at  all  frequencies  between  W  and  1/(2T). 

The  autocorrelation  sequence  of  the  sampled  process  may  also  be 
expressed  in  terms  of  the  normalized  spectral  density  function: 
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^(nT)  -  J  s3(u) 
-1/2  ' 


i2irnu  . 
e  du 


s3(u) 


-( 


Ts(u/T) ,  |u|  <  TW 

0  ,  TW  <  |u|  <  1/2 
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IV.  ANALYSIS  OF  InP  PHOTOCONDUCTOR  SWITCHES 
FOR  SIGNAL  SAMPLING 


by  T.O.  Jones  and  W.H.  McKnight 

This  analysis  is  intended  to  offer  a  first  consideration  of  signal 
sample  accuracy  (dynamic  range)  limitations  arising  from  the  physical 
properties  of  a  photoconductor  sample -/track- and-hold  switch  and  the  other 
circuit  elements  in  the  sample-/track-and-hold  circuit.  For  the  purpose  of 
this  analysis,  we  will  assume  that  a  typical  sample -/track- and-hold  circuit 
consists  of  four  basic  elements:  an  input  preamplifier,  the  photoconductor 
switch,  a  hold  capacitor,  and  an  output  amplifier.  Figure  10  illustrates 
the  corresponding  circuit  schematic  where 

Rq  ■  output  impedence  of  analog  preamplifier 

V(t)  ■  signal  source  to  be  sampled  (e.g.,  Vq  sin  wt) 

Rg(on)  *  resistance  (steady  state)  of  photoconductor  switch  when 
illuminated  (switch  closed) 

Rg(off)  ■  resistance  (steady  state)  of  photoconductor  switch  attribut¬ 
able  to  thermally  generated  carriers  (unilluminated  or  open- 
switch  state) 

Cg  =  capacitance  of  photocondutor  switch 

CH  ■  capacitance  of  hold  capacitor 

Ra  ■  input  resistance  of  sample  output  amplifier. 


Figure  10.  Sample  and  hold  schematic  circuit. 


This  circuit  functions  under  a  square-wave  clock  (e.g.,  mode-locked 
ion  argon  or  Nd  YAG  laser)  cycling  R„(t),  the  switch  resistance,  between  a 

b  7 

high-resistance  dark-state,  R  (off),  typically  >  10  ohms,  and  a  low- 

5  2 

resistance  illuminated  state,  R  (on) ,  typically  <  10  ohms.  Errors  in 
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accuracy  (dynamic  range  limitations)  occur  in  the  signal  data  samples 
produced  for  the  A/D  converter.  For  our  purposes,  we  will  consider  that  the 
(preamplified)  signal  voltage  ranges  from  -V  to  +Vq,  thereby  yielding  a 

quantization  step  size  of  2VQ/2n  -  VQ/2n  ^ .  Analysis  of  these  errors  is 

considered  in  several  categories:  (1)  tracking  error,  whereby  the  voltage 
across  C„  lags  V(t)  because  of  nonzero  resistance  R  +  R  (on) ;  (2)  discrete 

electronic  charge  error  resulting  from  statistical  fluctuations  of  charge 

flowing  onto  (or  from)  Cu  for  a  given  V(t)  time  profile  (with  R  closed); 

H  b 

(3)  pedestal/feedthrough  error  resulting  from  the  finite  turn-off  time  of 

the  photoconductor  switch  (finite  decay  time  of  optically  generated  charge 

carriers),  whereby  charge  will  flow  onto  (or  from)  C,,  after  the  laser  (clock 

fl 

pulse)  shuts  off,  including  feedthrough  error  from  charge  leakage  onto  (or 

from)  C  through  R  (off)  after  Re  reaches  a  maximum;  (4)  droop  error 
H  s  b 

resulting  from  discharging  through  RA  (when  Rg  is  open);  (5)  parasitic 

capacitance  error  resulting  from  the  nonzero  value  of  the  switch 
capacitance,  C  ;  (6)  thermal  noise  sources  in  the  resistive  and  capacitive 

b 

elements  of  the  circuit;  and  (7)  generation/recombination  noise  resulting 
from  the  generation  and  recombination  of  charge  carriers  in  the  switch. 

1.  Tracking  Error.  The  sample- /track- and-hold  cycle  starts  with  the 
laser  turning  the  photoconductor  switch  on  (illuminated  or  closed  switch)  at 
t  -  0  for  one-half  cycle  of  fg  (this  is  typical  in  A/D  converter  circuits), 

the  sample  frequency.  This  allows  the  hold  capacitor,  C  ,  to  track  the 

n 

signal  by  being  charged  through  R  (on)  and  R  .  The  current  through  C  is 

bO  b 

negligible  compared  to  the  current  through  R  (on)  and  the  turn-on  response 

b 

time  of  the  switch  is  negligible  (instantaneous)  compared  to  the  sample 
clock  period.  At  the  end  of  the  track  cycle,  the  voltage  uncertainty  or 
error,  AV^,  across  must  be  within  one-half  the  least  significant  bit 

(LSB)  of  the  "true"  voltage  (the  fractional  error  in  charge  on  C„  must  be 

less  than  one-half  the  LSB): 


VV2>  ■  6  -  2n+1 

where  vH(rs/2)  ■  true  signal  voltage  across  at  t  -  rg/2  (end  of  tracking 

period),  r  *  l/(f  ),  and  RC  ■  [R  +  R  (on)]C  .  One  can  solve  for  the  mini- 
b  b  O  b  n 

mum  resistance,  R  [where  In  2n+i  -  (n  +  l)ln  2  »  0. 7(n  +  1) ] . 


2.  Discrete  Electronic  Charge  Error.  For  the  discrete  electronic 
charge  error,  one  considers  that  the  number  of  electrons,  m,  that  flow 
through  the  switch  in  seconds  obeys  a  Poisson  distribution  with  an 

average  error  of  yin.  The  corresponding  voltage  uncertainty,  AV  must  be 

less  than  one -half  the  LSB: 


2C  V 
H  o 


■> 

max 


q  7m  2V 

o  max  o 

CH  -  2n+l 


where  is  the  electron  charge. 

Note:  During  a  span  of  r  seconds  (from  the  end  of  one  hold  period  to  the 

end  of  the  next),  the  greatest  change  in  signal  voltage,  V(t),  giving  the 
largest  flow  of  charge,  onto  occurs  for  the  highest  frequency 

component,  f  /2 ,  present  in  the  sampled  signal  by  an  amount  2V  (if  the 

O 

signal  contains  V  sinirf  t) .  Thus 

O  j 


This  gives  a  minimum  value  for  the  hold  capacitance. 


3.  Pedestal/Feedthrough  Error.  For  analysis  of  the  pedestal  error, 
one  considers  that  the  conductivity  of  the  switch  decays  exponentially  when 
the  light  is  turned  off  because  of  the  decay  of  the  optically  generated 
charge  carriers.  This  finite  decay  time  and  the  noninfinite  value  of 
Rg(off)  permit  charge  to  leak  from  or  onto  during  the  open- switch 

(unilluminated)  period  (from  t  -  r  / 2  to  t  -  r  ). 

J  J 

The  resistance  through  which  charge  leaks  from  (or  onto)  can  be 
modeled  according  to  the  following  schematic: 


where  Rc(t)  -  R  (on)  for  0  <  t  <  r  /2 
j  j  s 


2  r  t/r 

_  ,  .  o  o 

R  (on)e  e 


for  r  / 2  <  t  <  r( 

l)  k 


repeated 
every  sample 
cycle 


and  r  =  average  lifetime  of  optically  generated  charge  carriers  where  r  » 

O 

r  .  The  charge  leakage,  AQ  ,  during  the  time  the  switch  is  open  (r  /2  <  t  < 
O  H 

r  )  is  the  sum  of  charge  leakage  through  each  path: 


p  V  sinjrf  tdt  p  v  simri-cot  p 

H  ”  J  R  +  R(t)  +  J  R.(off)  <  VfS  J 


V  sirorf.tdt 
o _ o 


o  R  +  R_(on)e 
o  o 


*fSRS(°ff) 


3  2 

V  *  f  r  V 

o  So  _ o 

AQH  12R  +  wf _R_ (of f ) 

o  i>  a 

where  we  have  set  R  -  R  (on) ,  and  used  V  simrf  t  <  V  *f  t  and  a  simplifying 
time  shift. 

The  corresponding  fractional  voltage  change  is 

3  2 

AQh  *  fsro 

2VoCH<24RoCH  +2.CHRs(off) 

where  each  component  could  be  set  to  less  than  one-half  the  LSB. 

3  2 

n  f  r  ,  3R  C 

Pedestal  error:  __  0  <  - r  ->  r  <  ■  ,  ° - r 

3EoCh  -  2n-l  °  -  ,3  n-2 


t.  .  2nr 

Feedthrough  error:  „w„f£)  <  ^  ->  Rs<of£)  i  ^  . 

A.  Droop  Error.  At  the  beginning  of  the  hold  period  (at  t  -  r_/2), 
the  maximum  voltage  that  can  occur  across  is  V  ,  which  falls  according  to 


V:'.’  V  V  V 
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V2-* 

VH(t)  “V  A  H  ’  (Ts/2  ~  C  “  rS)  • 


The  voltage  across  C  at  t  -  r  seconds  is  thus 

H  b 


2rach 

VV  -  V 

so  that  the  maximum  voltage  error  would  be 


AVH  -  VV2>  -  W  “  vo  •  V 


2rach 


V  II  -  e 
ov 


2rach 


The  maximum  fractional  error  in  voltage  across  C  after  r  / 2  seconds 
would  be 


-  1  1  .  c2RaCh  <  -3— 

2Vo  2  2n+1 


5.  Parasitic  Capacitance  Error.  When  CH  has  been  charged  by  an 

amount  (representing  a  signal  value  V^)  and  the  sample  switch  opens,  some 

charge,  Qg,  on  will  flow  onto  (creating  a  voltage  Vg)  because  of  the 

changing  signal  voltage,  V(t),  which  will  thereby  change  the  voltage  across 
Cjj  by  an  amount  AV^.  Since  a  signal  of  V(t)  -  VQsin»rfgt  can  change  by  a 

maximum  of  J2V  during  r  / 2  seconds  (the  hold  period),  we  have  V(r  / 2)  = 

-V  /Jl  and  V(r  )  -  V  /Ji.  Thus  (for  perfect  tracking), 

O  bo 


VV2>  -  V<V2>  -  -±  -  5?  _>  Q"  ' 


V<V  -  VV  *  W  - 


Qe  V 
'  cs 


->  Q  -  -72V  r  I  r" 
S  o  Cg  +  CH 
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Qu  Qs  Vo  _ 

AV„  -  Vu(rc)  -  VH(rs/2)  -  7T  -  7T  +  ~  J2V, 


'H 


H  S 


CH  CH  '  f2 


CS  +  CH 


2V  “  J2 
o  y 


CS  +  CH 


<  r 
~  2n+l 


CS^ 


^CH 

2n+1  -  72 


6.  Thermal  Noise.  Thermal  noise  in  a  resistive  circuit  element  is 
typically  represented  as: 

<vth>  - 4kTRi£ 


where  k  -  Boltzmann's  constant;  T  ■  absolute  temperature;  R  ■  resistance  and 
Af  ■  noise  bandwidth.  This  open- form  expression  might  suggest  that  the 
noise  voltage  becomes  infinite  as  bandwidth  or  resistance  increases  without 
bound.  The  truth  is  that  this  formula  is  only  valid  within  some  limited 
bandwidth  beyond  which  the  noise  becomes  negligibly  small,  because  the 
parasitic  capacitance  associated  with  a  resistive  element  serves  as  a 
feedback  regulator  (shunt  capacitor)  for  noise  voltage.  For  an  infinitely 
wide  bandwidth,  the  thermal  noise  voltage  would  be  [8] 


kj 

C 


where  C  ■  parasitic  capacitance  associated  with  the  resistance.  Typical 
values  for  our  case  would  be  T  -  300“K,  C  -  lpf  so  that 

2  x  10  **  volts. 

Since  detection  of  a  microvolt  relative  to  a  volt  represents  a  dynamic  range 
of  20  bits,  this  noise  source  should  present  no  data  sample  accuracy 
limitations . 

7.  Generation/Recombination  (GR)  Noise.  This  is  a  noise  source  that 
arises  from  the  generation  and  recombination  of  charge  carriers  for  which 
numerical  values  of  the  spectral  density  function  of  the  current  fluctua¬ 
tions  typically  require  a  value  for  the  average  number  of  carriers  generated 
per  second.  Also  needed  are  values  for  the  average  lifetime  of  a  photo- 
excited  carrier,  the  charge  carrier  drift  mobility,  the  photoconductor 
switch  gap  length,  and  an  effective  voltage  across  the  gap.  Since  these 
parameters  can  typically  vary  by  a  considerable  extent,  further  considera¬ 
tion  will  be  left  to  the  future. 

It  should  be  mentioned  in  passing  that  the  concept  of  shot  noise 
(resulting  from  the  discrete  nature  of  electronic  charge)  typically  depends 
on  many  of  these  same  parameters  and  has  been  shown  [8]  to  be  related  to  GR 
noise  according  to 
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where  tq  ■  average  carrier  lifetime 

•  average  carrier  transit  time  across  gap. 

Typical  parametric  values  for  our  application  show  that 

*  200,  or  that  GR  noise  is  the  predominant  contribution  between 

these  two. 

The  foregoing  analysis  is  intended  to  be  a  first-cut  look  at  con¬ 
straints  and  limitations  arising  from  these  various  circuit/switch  param¬ 
eters  and  phenomena.  A  more  basic  first-principles  software  model  of  this 
situation,  involving  solutions  of  simultaneous  differential  equations  for 
charge  and  charge  flow  in  various  parts  of  the  circuit,  is  being  developed 
for  a  source  follower  configuration  in  which  the  FET  gate  capacitance  serves 
as  the  hold  capacitor.  An  empirical  effort  addressing  these  same  issues,  as 
well  as  the  impact  of  phase  noise  arising  from  clock  pulse  jitter  and  pulse - 
to-pulse  timing  errors,  is  also  under  way. 


V.  CURRENT  EFFORT  (FT  85-86) 

Present  development  is  proceeding  on  several  fronts: 

1.  The  linear  predictive  coding  technique,  as  described  in  [7]  and 
[8],  is  being  implemented  in  software  for  testing  and  analysis  for  suit¬ 
ability  in  this  application.  Critical  issues  include  the  impact  of  limited 
past  signal  samples  and  nonstationary  signal  statistics  on  prediction 
accuracy  and  stability. 

2.  The  impact  of  phase  noise,  arising  from  timing  pulse  jitter,  and 
mistiming  on  digitization  accuracy  is  being  empirically  investigated  by 
constructing  a  suitable  audio  frequency  (less  than  10-kHz)  sample-and-hold/ 
analog-to-digital  converter  system  and  imposing  amplitude-controlled,  band- 
limited  white  noise  or  appropriate  sinusoid(s)  on  the  clock  pulses.  This 
circuitry  would  utilize  conventional  electronic  sample-and-hold  switching 
with  16-bit  A/D  converters,  and  would  create  measurable  random  or  systematic 
errors  in  sample  pulse  timing  and  pulse  width  (duration).  These  effects  can 
be  measured  directly  as  to  their  impact  on  dynamic  range  limitations  of  the 
signal  being  digitized,  and  accordingly  scaled  to  actual  sample  frequencies 
on  the  order  of  100  MHz. 

3.  Development  of  circuitry  such  as  that  depicted  in  figure  11,  with 
a  state-of-the-art  commercially  available  A/D  converter  (offering  6  bits  of 
dynamic  range)  to  operate  at  100  MHz  is  being  pursued.  This  will  employ  an 
InP  picosecond  photoconductor  switch  fabricated  for  this  purpose  and  a  suit¬ 
able  diode  laser.  Thus  one  can  empirically  address  the  analysis  previously 


described  in  section  IV,  as  well  as  the  technique  employed  in  paragraph  2 
above,  at  the  actual  sample  rates  ultimately  to  be  used.  This  is  being 
carried  out. 


-100-MHz 
SAMPLING  RATE 


Figure  1 1 .  Track  and  hold  circuit. 


4.  InP  photoconductor  switches,  which  were  obtained  from  MIT  Lincoln 
Laboratories  in  November  1983,  have  been  found  to  exhibit  anomalous  behavior 
similar  to  that  observed  and  reported  in  the  fall  of  1982  [8]  (i.e.,  off- 
state  or  open-switch  resistance  of  100  ohms  or  less  after  initial 
illumination).  Efforts  are  currently  being  made  to  understand  and  deal  with 
these  observed  phenomena  as  well  as  the  relatively  low  observed  durability 
of  switch  lead  attachment.  These  efforts  include  switch  fabrication, 
characterization,  and  optimization  at  the  NOSC  microwave  integrated  circuit 
facility . 

5.  A  version  of  the  sample/track-and-hold  circuits  depicted  in 
figures  10  and  11  is  being  modeled  in  software  whereby  time  solutions  will 
be  obtained  to  the  simultaneous  differential  equations  governing  charge  and 
charge  flow  in  this  situation.  Thus  a  more  nearly  exact  time  response  can 
be  obtained  as  an  improvement  to  the  analysis  in  the  foregoing  Section  IV. 
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